Supersonic baryon-CDM velocities and CMB B-mode polarization 
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It has recently been shown that supersonic relative velocities between dark matter and baryonic 
matter can have a significant eflect on formation of the first structures in the universe. If this effect 
is still non-negligible during the epoch of hydrogen reionization, it generates large-scale anisotropy 
in the free electron density, which gives rise to a CMB B-mode. We compute the B-mode power 
spectrum and find a characteristic shape with acoustic peaks at £ m 200, 400, . . .. The amplitude of 
this signal is a free parameter which is related to the dependence of the ionization fraction on the 
relative baryon-CDM velocity during the epoch of reionization. However, we find that the B-mode 
signal is undetectably small for currently favored reionization models in which hydrogen is reionized 
promptly at 2 ~ 10, although constraints on this signal by future experiments may help constrain 
models in which partial reionization occurs at higher rcdshift, e.g. by accretion onto primordial black 
holes. 

PACS numbers: 



I. INTRODUCTION 

It was recently realized [1] that supersonic relative ve- 
locities between baryons and cold dark matter can have 
a significant effect on formation of the first structures in 
the universe. The physics can be described intuitively as 
follows. 

The baryon-CDM relative velocity W6c(x) is coherent 
on scales smaller than a few comoving Mpc (Fig. 1). Let 
us imagine dividing the universe into subregions of this 
size, with a constant value of vtc in each subregion. As 
shown in [1], the baryon-CDM flow will suppress growth 
of perturbations whose comoving wavenumber k is larger 
than the freestreaming scale kps = aH/vbc- Equiva- 
lently, a mode is suppressed if its wavelength is less than 
the distance that the baryon-CDM flow travels in a Hub- 
ble time. 

Perturbations whose wavenumber k is larger than the 
Jeans scale kj = aH/cg, where Cg is the baryon sound 
speed, are additionally suppressed by the usual Jeans 
mechanism (pressure support). For large fc, Jeans sup- 
pression dominates (by an extra factor of k in the equa- 
tions of motion) over the suppression. Therefore, 
the Vbc suppression is important on a window of scales 
given by kps < k < kj, illustrated in Fig. 2. The size 
(kj/kps) of this window is equal to A4 — Vbc/cs, the 
Mach number of the flow, so the existence of a range of 
scales where the Vhc suppression dominates is equivalent 
to a supersonic flow. For z < 100, a typical Mach num- 
ber is M = {vl^^^"^ /cs ~ 1.5, so there exists a range of 
scales where Vbc suppression is an order unity effect. 

The first stars in the universe are predicted to form 
at redshift z ^ 30, in dark matter halos with typical 
mass M ^ 10^ Mq. The abundance of such halos de- 
pends sensitively on the amplitude of perturbations on 
a broad range of scales near characteristic wavenumber 
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FIG. 1: Power spectrum of the baryon-CDM relative velocity 
field at z = 10.5. The power spectrum drops sharply 
for fc > 1 /i Mpc~^, indicating that the baryon-CDM flow is 
coherent on scales larger than a megaparsec. The shape of the 
power spectrum is roughly independent of z (for z < 100), but 
its amplitude varies as A^i,^{k) = Pvbc{k) /i-K^ oc (1 + 2)^, 
with RMS velocity (wL)^''^ = 0.36 Km/s at 2 = 10.5. 



k = (47rp„/3M)i/3 - 80ft. Mpc"^ This range of scales is 
broad enough that there is some overlap with the range 
kps <k <jkj where v^c suppression is important, and so 
the abundance of the first stars is sensitive to the local 
value of Vbc- For larger halo masses, the overlap of scales 
decreases rapidly, and vi,c suppression is less important. 

If we spatially average over a large region of the uni- 
verse, then the overall effect of the baryon-CDM flow is to 
delay the formation of the first bound structures. What 
is more distinctive is that the effect is inhomogeneous: 
the halo abundance will depend on the local value of Vhc- 
Since the power spectrum of the v^c field has a differ- 
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observables. For example, a recent paper [2] considered 
the possibility that the Vbc effect survives long enough to 
affect galaxy abundance at redshifts z <b. Here, we will 
study the analogous possibility that the v^c effect survives 
long enough to affect the hydrogen ionization fraction 
during the epoch of reionization (z ^ 10), and compute 
the CMB B-mode which is generated. In particular, we 
would like to understand whether large-scale B-modes 
sourced by the Vhc effect can significantly contaminate 
the gravity wave signal at ^ ~ 50, since this will be a 
primary scientific target for polarization experiments in 
the near future. 



II. B-MODE CALCULATION 



FIG. 2: Redshift evolution of length scales relevant for the 
Vbc effect. As described in the text, the baryon-CDM flow 
suppresses structure formation on a range of scales given by 
kps < kj. This makes the abundance of halos of mass 
< 10® Mq sensitive to the local value of vi,c, but the effect 
decreases quickly with increasing halo mass. 



ent shape than the usual contributions to the halo power 
spectrum (i.e. a Poisson term and a term proportional 
to the matter power spectrum), this makes the v^c effect 
qualitatively new, e.g. it was shown in [2] that the Vbc 
effect can shift the apparent scale of the baryon acoustic 
peak in the halo two-point function. 

At redshift z ^ 30, the abundance of halos (or stars) is 
significantly modulated by the local value of v^c- What is 
less clear (and a critical question for the phenomenology) 
is whether the effect "survives" to lower redshift: if we 
inspect the universe at z = 10, can we still tell the differ- 
ence between a low-Wf,c and a high-Ufcc region? The gravi- 
tational suppression mechanism described above suggests 
that the abundance of the largest halos is relatively unaf- 
fected by Vhc'- the scales which are suppressed by Vbc are 
smaller than the scales which are relevant for spherical 
collapse of these halos (Fig. 2). However, the number 
(or total luminosity) of luminous objects such as galaxies 
could depend on vtc, due to a variety of possible feedback 
effects. An example of positive feedback would be metal 
enrichment of the intergalactic medium: in a large-Whc 
region, formation of the first stars is suppressed, leading 
to lower IGM metallicity, which further suppresses star 
formation at later times since it will be harder for gas to 
cool via atomic lines. In the opposite direction, an exam- 
ple of negative feedback would be heating of the IGM by 
starlight: in a large-Ubc region with fewer stars at z ~ 30, 
the IGM will be colder than average, which makes it eas- 
ier for clouds of gas to overcome pressure support and 
collapse gravitationally, leading to increased star forma- 
tion at later times. 

It is also interesting to take an agnostic approach to the 
question of whether the V{,c effect survives to low redshift 
(i.e. treating the size of the effect as a free parameter), 
and study the potential impact on various cosmological 



We assume that the mean ionization fraction Xe{z) at 
redshift z depends weakly on the local value of v^c- By 
rotation invariance, Xe can only depend on v'^^; it is con- 
venient to change variables by defining the field 



""be 



(r,z) 



"be 



(1) 



which is normalized to have zero mean, and is (to an ex- 
cellent approximation) independent of redshift as implied 
by the notation. 

We will assume that the Vbc effect is small at z ^ 10, 
so that the ionization fraction Xe{z) can be linearized in 
■q: 



Xe(r, z) = Xeiz) (1 + b^e{z)'n{^)) 



(2) 



where we have introduced a free parameter b^eiz) which 
parameterizes the size of the Vbc effect at redshift z. In- 
tuitively, bxe{z) is the mean fractional change in Xe{z) 
produced by a la fluctuation in the local value of Vbc- 
The CMB optical depth in direction n is given by 



r(ri) = ne^oCT J dz 



(1 + ^)^ 
Hiz) 



Xe{xiz)n,z) (3) 



where n^.o is the comoving electron density at redshift 
zero, ax is the Thomson scattering cross section, and 
x(z) is the comoving distance to redshift z in the assumed 
flat cosmology. In the model (2), the angular power spec- 
trum of the r field is given (in the Limber approximation) 
by 
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Note that we are only considering new contributions to 
CJ'^ due to variations in the ionization fraction due to 
the Vbc effect; the full power spectrum contains addi- 
tional terms (e.g. a "one-bubble" term arising from Pois- 
son statistics of individual HII regions) not considered 
here [3]. 
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Linear CMB polarization is generated by Thomson 
scattering of CMB photons by free electrons. If reioniza- 
tion is spatially homogeneous, then the combination of 
rotation and parity invariance implies that only E-mode 
polarization is generated. The CMB B-mode power spec- 
trum generated by inhomogeneous reionization is a sum 
of scattering and screening terms [4-6]; these correspond 
respectively to generation of new linear polarization via 
Thomson scattering by free electrons, and screening of 
the primary E-mode by inhomogeneities in the optical 
depth. 
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where t^s ~ 0.09 is the effective optical depth to reion- 
ization and Qrms = 17.7 /xK is the RMS temperature 
quadrupole during matter domination. 

For calculating the B-mode power spectrum at the 
wlO% level, we can make the approximation that the 
comoving distance xl^^) is constant over the range of red- 
shifts which contribute to the integral in (4). We can 
then factor CY as the product of a parameter a which 
depends on reionization history and the bias parameter 
hxe^z), and a template shape in I which depends only on 
cosmological parameters: 
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To complete the calculation of the B-mode power spec- 
trum, we need to compute the power spectrum Pr,{k). 
The two-point correlation function of the vector field vl^ 
is given by: 

hi h.3 

Kc(k)*<(kO) = -^P„6c(fc)(27r)3^3(k - k') (11) 

where i,j are spatial indices and Pvbc{k) is the relative 
velocity power spectrum, which can be computed using 
CAMB [7]. From the definition (1) of 77, the power spec- 
trum Prj{k) is a four-point function in Vbc which can be 



computed straightforwardly using (11) and Wick's theo- 
rem. The result can be presented cither in Fourier space, 



P,(fc) = 
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or equivalently in position space as 

P^{k) = 47r / drr^joikr) {6Mrf + ^Mr)^) (13) 



where the correlation functions tpoir), ''p2{r) are defined 
by: 

1 f h"^ dh 

Mr) = -^Pvbc{k)jo{kr) (14) 
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In Fig. 3, we show the B-mode power spectrum for a 
fiducial reionization model defined as follows. The reion- 
ization history Xe(z) is given by the CAMB parameteri- 
zation, with reionization width Az = 1.5, central redshift 
Zrei = 10.5, and total optical depth r = 0.09. To get a 
sense for the order of magnitude of a, if hxe{z) is con- 
stant in z, then a = 0.00866^^. In [8], the value of 6^^ at 
2; ~ 10 is predicted to be ~ 0.01 (this value is consistent 
with the first attempts to include the v^c effect in sim- 
ulations [9-12]), so we will take a = 8.6 x 10~^ as our 
fiducial value. 

For this fiducial reionization model, we find that the to- 
tal B-mode power spectrum generated by baryon-CDM 
relative velocities is small (roughly 0.14 /uK-arcmin on 
large scales). The scattering power spectrum Cf""^ has 
acoustic peaks at values of £ which are multiples of 
^* = TTXrei/sBAO) where Xrci ~ 9800 Mpc is the angu- 
lar diameter distance to reionization and sbao ~ 150 
Mpc is the baryon acoustic scale. 

There is one additional c^ffcict which has not been in- 
cluded in our B-mode calculation: high-£ damping due 
to freestreaming. Our model (2) allows the mean ion- 
ization fraction x^, at position r to depend on the veloc- 
ity f6c(r) at the same point r; in reality it will also de- 
pend on nearby values of Vbc since the mean free path of 
reionizing photons is nonzero. The qualitative effect will 
be to suppress the power spectra C^""" , C^"" on scales 
£ > (xroi/A) « 1000, where A « 10 Mpc is the char- 
acteristic size of HII regions during reionization, but a 
precise calculation of the damping tail will depend on 
model-dependent details. 



III. OBSERVATIONAL PROSPECTS AND 
DISCUSSION 

To forecast detectability of the CMB B-modes pro- 
duced by the Vbc effect, we consider several observational 
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FIG. 3: B-mode power spectra ""^ and generated 



by baryon-CDM relative velocities, with the lensing B-mode 
power spectrum shown for comparison. For I > 1000, the 

B-mode power spectra Cf""^ , Cf'" will be additionally sup- 
pressed (beyond what is shown here) by photon freestreaming 
during the epoch of reionization. For I < 1000, these power 
spectra have shapes which are roughly independent of the 
reionization model, but the overall amplitude a is a free pa- 
rameter. In the figure, we used a — 8.6 x 10~^, corresponding 
to the fiducial reionization model described in the text. 



scenarios. The la error on a is given by 



cr(a) 



-1/2 



(16) 



where is the lensing B-mode power spectrum and 

Nf^ is the noise power spectrum of the experiment. 

For a lensing-limited experiment (i.e. Nf^ <^ C^'""'') 
with /sky — 0.7 and iamx = 1000, we find a{a) — 
8.6 X 10"''. This forecast effectively treats the lensing B- 
mode as an extra source of noise and could potentially be 
improved by incorporating delensing techniques (e.g. [IS- 
IS]) which statistically separate the lensing B-mode from 
other components. Using the forecasting methodology 
from [16], we find that the statistical error can be im- 
proved to a{a) = 7.5 x 10"'*, 2.7 x 10"'', or 1.3 x 10"^, 
assuming a Gaussian beam with 6'fwhm — 3' and po- 
larization noise level Ap = 3, 1, or 0.25 /iK-arcmin re- 
spectively. In all cases, the statistical error cr{a) is larger 
than the value aM = 8.6 x 10^^, indicating that no sta- 
tistically significant detection of the Vhc signal is possible 
in our fiducial reionization model. 

The basic reason that the Wf,c B-mode is small in our 
fiducial model is that we have assumed that reionization 



is dominated by stars at z < 10, and that the abundance 
of these sources depends weakly on the local value of Vbc 
(i.e. bxe is small at these redshifts). A larger B-mode 
may be generated if bxe can become large. For example, 
if positive feedback mechanisms (such as metal enrich- 
ment of the IGM) help the effect survive to z ^ 10, 
then may be enhanced, leading to a larger signal. 
Models of reionization in which high redshifts contribute 
non-negligible optical depth should also lead to increased 
bxe- For example, scenarios where X-ray emission from 
primordial black holes ionize the universe by a few per- 
cent at redshift ^ 100 have been considered e.g. in [17]. 
Although a quantitative study is beyond the scope of 
this paper, the bias parameter bxe could be as large as 
order unity in these models, since the Hoyle-Bondi rate 
for accretion of gas onto a black hole is proportional to 
(ff,c + Cs)^^/^, and and are of the same order of 
magnitude (so that statistical fiuctuations in Vhc lead to 
order- unity variations in the accretion rate). 

In conclusion, although optical depth anisotropy 
sourced by baryon-CDM relative velocities generates a 
CMB B-mode in principle, we find that the amplitude 
is likely to be unobservably small, if currently favored 
reionization models are correct. On the other hand, no 
new data analysis is required to fit for the amplitude a 
of the Vhc B-mode (one simply fits a multiple of a fixed 
template shape to the estimated power spectrum) , so es- 
timating a from future B-mode measurements may be a 
useful test of the assumptions underlying these models. 
In any case, since the acoustic features at i = 200, 400, . . . 
are on different scales than the ones which are relevant 
for gravity waves (£ < 10 and £ ~ 50), the B-mode power 
spectrum calculated here cannot be a contaminant for the 
gravity wave signal, even in a reionization model where 
its amplitude is large. 
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